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Summary. We have developed a quick and accurate numerical tool to compute 
gas cooling whichever its chemical composition. 

1 Introduction 

Metal abundances in galaxies vary from place to place and with time. These 
variations reflect the star formation history within galaxies because the chem- 
ical enrichment results from synthesis of heavy elements by successive gen- 
erations of stars. Gas radiative cooling is very sensitive to the ISM chemical 
composition. Cooling influences both the gas dynamics and the star forma- 
tion rate. Therefore we need to follow self-consistently dynamical and chemi- 
cal evolution. We have developed a method to compute in a very short CPU 
time accurate cooling functions dependent on the exact abundances of el- 
ements. It will be used in high resolution chemodynamical simulations of 
galaxies. 

2 Description of the method 

It is well known that cooling rates depend on gas chemical composition. A gas 
with solar composition has a cooling rate greater than without metals by more 
than an order of magnitude for some temperatures (Fig. 1). Cooling rates are 
generally interpolated between rates computed for different metallicities with 
solar abundance ratios (see [1] for instance) or with solar abundance ratios 
and enhanced abundances for some elements (see [2] for instance). The gas 
abundances are obviously not always solar because the gas ejected from stars 
by winds and SNe has a non-solar composition. Hence we need to take into 
account the real abundances of each element in the gas to build more realistic 
cooling functions whichever the chemical composition. However, computing 
cooling rates with Mappings III for each gas particle at each timestep of a 
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Fig. 1. Cooling rates as a function of temperature for different metallicities com- 
puted with Mappings III. The abundances used in calculations are taken from [2] 



chemodynamical simulation is a very CPU time consuming task. We thus 
have developped a recipe to reconstruct the cooling curve on the fly. 

We consider the ISM as an optically thin gas in collisional ionization equi- 
librium. Calculations are performed in the temperature range 10 4 K to 10 8 K 
using Mappings III, successor of Mappings II whose cooling computations 
details are described in [2]. All the abundances of elements are relative to 
hydrogen. The hydrogen density is 1.0cm -3 , independent of temperature. 
We drop the temperature dependence of cooling rates (A(T)) for the sake of 
legibility. 

Many cooling processes depend on the density of electrons and on the 
abundances of elements. A simple hypothesis is to suppose that cooling rates 
are proportional to these two quantities. Thus, we can try to reconstruct the 
total cooling curve A to t with a linear combination of individual curves for 
each element. We first need to compute these individual curves. We obtain 
Ah with a purely hydrogen gas. For all the heavier elements used for calcula- 
tions in Mappings III (He, C, N, O, Ne, Na, Mg, Al, Si, S, CI, Ar, Ca, Fe, Ni), 
we have calculated cooling rates for mixtures of hydrogen and one element Xi 
at its solar abundance Ax u h- Individual cooling curves Axi have then been 
obtained substracting Ah from each An,Xi- Ax i = An.Xi — Ah- For a given 
chemical composition, we then compute the quantity «x, for all the elements 
previously listed, where aj ; is the abundance of element Xi normalized by 
its solar abundance: aj, = nxJnx iQ - For helium, our hypothesis is too sim- 
ple. Its cooling rate cannot be considered proportional to its abundance with 
a good approximation as there are too strong non-linearities. We choose to 
linearly interpolate in a grid of precomputed curves for different abundances 
to obtain the helium cooling curve Ane,a H ^ ■ As hydrogen and helium are the 
two most abundant elements, we assume that they are the only sources of 
electrons. We indeed consider that the contribution of the metals are negligi- 
ble because of their much lower abundances in the ISM. By construction, our 
individual cooling functions for metals take into account the electrons of the 
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Fig. 2. (left panel) Cooling rates for several elements at solar abundances. The 
curve for other elements is the sum of the contributions of N, Ne, Na, Mg, Al, Si, S, 
CI, Ar, Ca and Ni. The solid broad line Atot is the cooling curve for solar abundances 
reconstructed using our recipe, (right panel) Percentage of relative error made 
when reconstructing cooling functions with our recipe for different compositions 

clement and hydrogen. Thus we only need to add the contribution of electrons 
provided by the ionization of helium ne (He). We add the terms for hydro- 
gen, helium and metals to build Atot so that we finally obtain the following 
approximation: A to t ~ A H + A H e,a Ha + (1 + ne (He)) Y^metais a x,A Xi 

To check the accuracy of our tool, we have compared cooling curves com- 
puted with Mappings III and the ones reconstructed with our recipe for 
various chemical compositions. We choose a solar composition to test a well- 
known chemical composition. Carbon, oxygen and iron are three of the most 
important coolants (Fig. 2, left panel). Furthermore, they are three of the 
most abundant elements in stellar winds and supernovae ejecta. Thus, we 
increase by a factor of 10 the abundance of these three important coolants. 
This factor is somewhat arbitrary and was only chosen to test the accuracy 
of our recipe in the case of a very high-enhanced gas. In all cases, the relative 
errors remain below a few percent (Fig. 2, right panel). The errors due to the 
computation of cooling rates are thus comparable to the others coming from, 
for instance, the hydrodynamical scheme. 
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